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Stripe orders in the extended Hubbard model 
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We study stripe orders of charge and spin density waves in the extended Hubbard model with 
the nearest-neighbor Coulomb repulsion V within the mean field approximation. We obtain 
V vs. T(temperature) phase diagram for the on-site Coulomb interaction U/t = 8.0 and the 
filling n = 0.8, here t is a nearest-neighbor transfer energy. Our result shows that the diagonal 
stripe spin density wave state (SDW) is stable for small V, but for large V the most stable state 
changes to a charge density wave-antiferromagnetic (CDW-AF) state. Especially we find at low 
temperature and for a certain range of value of V, a vertical stripe CDW-AF state becomes 
stable. 
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§1. Introduction 

Much attention has been focused on the phe- 
nomena of charge and spin orderings in transition 
metal oxides including high temperature superconduc- 
tiong cuprates since the discovery of the 1/8 prob- 
lem of La 2 - x Ba a ,Cu04. 1 ) La2- x Sr x Ni04 2 ~ 15 ' and 
I^NiCu+s 16-22 ) is the typical material of this phenom- 
ena. L^-zSr^NiC^ shows the stripe order of the spin 
and the charge and is insulating while x < O.5. 2 ' The di- 
rection of the stripe is diagonal to the crystal axis. 2 ) The 
relation between the hole concentration nh and the in- 
commensurability e is approximately linear. 5 ' The stripe 
order consists of the charge stripe and incommensurate 
spin density wave(SDW). The critical temperature of the 
charge stripe is higher than that of the incommensurate 
SDW. 5 ) 

In early theoretical mean-field study of the Hubbard 
model 23-27 ) and the d-p model, 28 ' the stable states are 
the insulating stripe or incommensurate spin density 
wave (SDW) states. In the hole density vs. on-site 
Coulomb interaction phase diagram at zero temperature, 
the diagonal stripe SDW state becomes stable for large 
on-site Coulomb intercation and small hole doping. 26 ' 
Temperature dependence of the stable states is studied 
within the mean- field approximation in Ref. 29. The 
result shows that as decreasing temperature, the stable 
state changes from the paramagnetic state to Neel state, 
the spiral SDW state, the stripe SDW state and the un- 
dulated stripe SDW state, successively. Especially, the 
charge and the spin stripe orders appear simultaneously. 
But this result is inconsistent with the experimental fact. 

Also such charged stripe may become unstable by the 
long range Coulomb interaction, although there is the 
distortion of the ionic lattice. Therefore we study the ex- 
tended Hubbard model with nearest-neighbor Coulomb 
interaction in order to consider the stability of the stripe 
state against the long range Coulomb interaction. We 



also consider the temperature dependence of the stable 
state within the mean-field approximation in order to 
search for the possibility of the charged stripe without 
the static spin ordering at finite temperature. 

This paper is arranged as follows. In section 2, we 
show our model and our method of the mean-field ap- 
proximation. In section 3, we show our numerical result, 
especially the phase diagram and we examine the stable 
state using the band calculation. Section 4 is devoted to 
the discussion and conclusion. 

§2. Model 

We use following the extended Hubbard model, 
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(2.1) 



where are nearest neighbor sites, cj {ci a ) is an 

electron creation (annihilation) operator at site i with 
spin (j, h irT = CfoCie (er =1,4) is an electron number 
operator at site i with spin a and n = n-r + n^. The 
nearest-neighbor transfer energy is t and the on-site and 
the nearest-neighbor Coulomb interactions are U and V, 
respectively. 

We take following mean fields for each site; 
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(2.3) 



where ( ) means thermal average. Also we take following 
mean fields for each nearest-neighbor bond 



Wirj*' = (clc ja ,) . (2.4) 

Then our mean- field Hamiltonian becomes as follows, 
T-Cmf = -t ( c L c iv + H - c ) + ^Upivcl^Ci-v 

<i,j> i,cr 
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where 



p t = Pi] + Pa- 
We write n-th eigenstate of this hamiltonian as, 

|vf„) = Eoilo}, 



(2.5) 



(2.6) 



(2.7) 



with the eigenvalue E n . Then Schrodinger equation be- 
comes, 

n.n.(i) 

-t Y, ^ + u Pl - a w a ~uxfw_ a 

3 

n.n.(i) n.n.(i) 

+V E PiC - V E ^V*^ = S„^,(2.8) 

where n.n.(i) means nearest neighbor sites of i-th site 
and X* means X*(Xi) when a is T(l)- We can rewrite 
mean fields using eigenstates {ip™ a } as follows; 

««r = X)l^| 2 /(^n-/i), (2-9) 

71 

x^E^a/O^-M), (2.10) 

= e c;^/ - M ) , (2.H) 

where / is the fermi distribution function and // is a 
chemical potential. Also we impose the electron number 
conservation, 



N e = E / (En - P) , 



(2.12) 



where N e is the total number of electrons. 

We solve eq. (2.8), then find out the chemical poten- 
tial p by eq. (2.12) and substitute the results to self- 
consistent equations, eqs. (2.9), (2.10) and (2.11). We 
continue this procedure until all of the mean fields con- 
verge. For each values of U, V and temperature T, we 
use several initial states and obtain several self-consistent 
solutions. Then we calculate the free energy of each state 
using following formula, 

f = pN e - fc B rE ln i 1 + e-n^-ri] - c/E^tm 



i m - 



t(TJ<7 



(2.13) 



Comparing free energies we determine the most stable 
configuration for each values of parameters. 



§3. Results 

We use a square lattice with 10 x 10 sites and im- 
pose the periodic boundary condition. The total electron 
number is fixed to 80, and this means that the filling is 
0.8. Also we fix the on-site Coulomb interaction U to 8, 
where we choose t = 1 as an energy unit. Form the re- 
sult of Ref. 26, the diagonal stripe state is stable at zero 
temperature and V = 0. We vary V from to 3 in units 
of 0.25 and T from 0.1 to 3.75. For each values of both 
of parameters V and T, we choose several initial states 
for the self-consistent procedure, and obtain several self- 
consistent solutions, which may be the global minimum 
or the local minimum of the free energy. 

3.1 Phase diagram 

In the following, we show the transition of most stable 
state as we vary the temperature for each value of V . 

When nearest neighbor Coulomb interaction V=0, de- 
creasing temperature T the most stable state changes 
from the paramagnetic state to the Neel state, then 
changes to a diagonal stripe SDW state. This result 
almost agrees with the previous result with U — 5, 29 - 1 
although in the previous result shows that the lowest 
temperature phase becomes undulated stripe SDW state 
which becomes stable at T < 0.1, the vertical stripe state 
is stable instead of the diagonal stripe state and the spiral 
SDW state becomes stable in a small temperature region 
between the Neel state and the diagonal stripe state. In 
Fig. 1, we show temperature dependence of free energies 
of several phases at V = 0. In this figure, we plot free en- 
ergy difference from the paramagnetic state. From this 




Fig. 1. Temperature dependence of free energies for several stable 
states when V = 0. Each free energy is measured from the free 
energy of paramagnetic state. 



figure we find that, the diagonal stripe state becomes 
more stable than the vertical stripe state when temper- 
ature is lower than a critical point where these stripe 
states become more stable than the Neel state, although 
the free energy of the vertical stripe state is lower than 
that of the diagonal stripe state in the region where the 
Neel state is most stable. We show the spin density, the 
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charge density of the diagonal stripe state at T = 0.1 and 
V = 0.0 in Fig. 2, the bond spin and the bond charge 
densities in Fig. 3. Here the charge density is pi and the 
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Fig. 3. Diagonal stripe state at V = and T = 0.1. a) Bond 
spin density, b) Bond charge density. The meaning of the arrow 
is same as that of Fig. 2. The size of the sphere shows the 
magnitude of the bond spin density. 



Fig. 2. Diagonal stripe state at V = and T = 0.1. a) Spin 
density, b) Charge density. The size (direction) of the arrow 
shows the magnitude (direction) of the moment. The sphere 
(tetrahedron) means the charge density is larger (smaller) than 
the average value of the charge density (0.8), respectively, and 
the size of the sphere (tetrahedron) shows the deviation from the 
average value. 



spin density is defined as follows, 

Pit 



), 



(3.1) 



where K and 3 mean the real and the imaginary part 
of a complex number, respectively. We define a complex 
bond charge density for a bond between nearest neighbor 
sites (ij) as, 



Pij 



lajcr- 



(3.2) 



Also we define the complex bond spin density as, 

Sij = 



(- 



w im + w, 



W, 



W, 



w, 



»Tj'T 



Wi 



iljl 



2i 



)(3.3) 



Note that ipi a is real, therefore pij, Sfj and SL are real 
and is pure imaginary. All of stable states that we 
have obtained do not have the imaginary Sfj component. 
So we plot real part of p^ and real part of Sy at the 
center of the bond (ij) in Fig. 3.. 

For the paramagnetic and the Neel states the bond 
charge density distribution is uniform and the bond spin 
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Fig. 4. Vertical stripe state at V = and T = 0.1. a) Spin 
density, b) Charge density. The meaning of symbols is same as 
that of Fig. 2. 
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Fig. 5. Vertical stripe state at V = and T = 0.1. a) Bond spin 
density, b) Bond charge density. The meaning of symbols is 
same as that of Fig. 3. 



Fig. 6. CDW-AF state at V = 1.0 and T = 0.5. a) Spin density, 
b) Charge density. The meaning of symbols is same as that of 
Fig. 2. 



densities vanish. In contrast to this, both of the bond 
spin and bond charge density distributions of the diago- 
nal stripe state have a diagonal stripe structure. These 
densities are large at the bonds that connect the site in 
a stripe and the site in an antiferromagnetic (AF) do- 
main. We also show the vertical stripe states in Fig. 4 
and Fig. 5. Similar to the diagonal stripe state, the bond 
spin and bond charge densities have stripe structure and 
are large at the boundary between of the stripe and the 
AF domain. 

For finite but small V (< 1.0), the temperature depen- 
dence of the most stable state is same as V = case. 
When V is increased, the structure of the diagonal stripe 
state is invariable, though the magnitudes of all of its 
densities become small. This means the diagonal stripe 
state is robust about small variation of the nearest neigh- 
bor repulsion V. Also the structure of the paramagnetic 
and the Neel states are invariant with respect to the vari- 
ation of V. 

Further increasing V up to 1.0, then a charge density 
wave-antiferromagnetic state (CDW-AF) becomes stable 
in a small temperature region which locates between the 
Neel state and the diagonal stripe state. We show the 
structure of the CDW-AF state in Fig. 6. Also in Fig. 7, 
we show the variation of free energies of several states 
with temperature. The charge density distribution of the 
CDW-AF state is same as the ordinary charge density 
wave state (CDW). However, the spin density does not 
vanish and is similar to that of the ferrimagnetic state; 
large moment appears at the site where the charge den- 
sity is large and opposite small moment appears at the 
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Fig. 7. Temperature dependence of free energies for several stable 
states when V = 1. Free energies are measured from the free 
energy of paramagnetic state. 



site where the charge density is small. The bond spin 
and bond charge densities are similar to the paramag- 
netic and the Neel states; the bond charge density is 
uniform and the bond spin density is very small. 

Further increasing V up to 1.5, the temperature de- 
pendence of the most stable state is almost same, but 
the temperature region where CDW-AF state is most 
stable, becomes large. When V = 1.75, a vertical stripe 
CDW-AF state becomes more stable than the diagonal 
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stripe state at low temperature (T = 0.1). Although, 
at T = 0.1 we obtain a vertical stripe SDW-CDW state 
which is a mixture of the vertical stripe SDW state and 
the CDW state, its free energy is slightly higher than 
that of the vertical stripe CDW-AF state. Also we have 
not obtained the diagonal type of stripe CDW-AF or 
stripe SDW-CDW state. In Fig. 8, we show temper- 
ature dependence of free energies of several states at 
V = 1.75. Also we show the spin and charge distri- 




T 




Fig. 8. Temperature dependence of free energies for several stable 
states when V = 1.75. Free energies are measured from the free 
energy of paramagnetic state. V-S-CDW-AF stands for vertical 
stripe CDW-AF and V-S-SDW-CDW stands for vertical stripe 
SDW-CDW. 



Fig. 9. Vertical stripe CDW-AF state at V = 1.75 and T = 0.1. 
a) Spin density, b) Charge density. The meaning of symbols is 
same as that of Fig. 2. 



butions of this new phase, the vertical stripe CDW-AF 
state, at V = 1.75 and T = 0.1 in Fig. 9, the bond spin 
and the bond charge distributions in Fig. 10. From these 
figures, we can consider this phase as a vertical striped 
previous CDW-AF state. From the point of view of the 
CDW order, this phase has vertical domains where the 
CDW order is large, which are separated by the stripes 
where the CDW order is small. However, from the point 
of view of the antiferromagnetic order, this phase has 
narrow domains of the ferrimagnetic phase, which is sep- 
arated by the broad stripes of a half ferromagnetic phase, 
where half of the sites have ferromagnetic order and the 
others have no spins. Similar to the stripe SDW state, 
the bond spin density is large at the bond between the 
ferrimagnetic domain and the half-ferromagnetic stripe. 
The bond charge density is large at the bond where bond 
spin density is large. However the bond at the both side 
of the site where magnetic moment is large in the half- 
ferromagnetic stripe, also has large bond charge density. 

In the parameter region V — 1.75 and 2.0, the temper- 
ature dependence of the stable state is almost same, but 
the temperature region where the vertical stripe CDW- 
AF state is stable becomes large at V — 2.0. 

At V = 2.25, the higher temperature phase changes 
from the Neel state to the CDW state. The spin and 
the bond spin densities vanish and as shown in Fig. 11, 
the bond charge distribution is uniform. This behav- 
ior can be expected, because the large nearest-neighbor 
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Fig. 10. Vertical stripe CDW-AF state at V = 1.75 and T = 0.1. 
a) Bond spin density, b) Bond charge density. The meaning 
ofsymbols is same as that of Fig. 3. 
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Fig. 11. CDW state at T = 1.5 and V = 2.25. a) Charge density, 
b) Bond charge density. The meaning of symbols is same as that 
of Figs. 2 and 3. 



Coulomb repulsion causes the alternating charge density 
distribution. The temperature dependence of free ener- 
gies of various states is shown in Fig. 12. 
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Fig. 13. Temperature dependence of free energies for several sta- 
ble states when V = 3.0. Free energies are measured from the 
free energy of paramagnetic state. 



than that of CDW based states. 

The total phase diagram the temperature T vs. the 
nearest-neighbor repulsion V is shown in Fig. 14. 
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Fig. 12. Temperature dependence of free energies for several sta- 
ble states when V = 2.25. Free energies are measured from the 
free energy of paramagnetic state. 



Fig. 14. Phase diagram in T vs. V plane. 



At V = 3.0, stripe phases are no longer stable and the 
most stable state changes from the CDW to the CDW- 
AF state as decreasing temperature. The temperature 
dependence of the free energies of various states is shown 
in Fig. 13. Free energies of stripe orders are much higher 



3.2 Energy band 

Here we discuss the electronic structure of several sta- 
ble states. In Fig. 15, we show the energy band of the 
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diagonal stripe states at T = 0.1 and V = 0.0 and 1.5. 
There are three types of the energy band. One is the oc- 
cupied valence band, second one is the midgap band and 
third one is the empty conduction band. The wave func- 
tion of the mid-gap state is localized at stripes. The site 
on the diagonal stripe does not have the nearest neighbor 
site on the stripe. Therefore the band is almost flat. The 
nearest neighbor Coulomb interaction V does not much 
affect the band structure, only the lower band become 
split around B-point(fc = (3tt/5, 2tt/5)). At V = 0.0 and 
at T-point(fc = (0,0)), the wave functions of the lower 
band of the valence bands have large amplitude on the 
edge sites of AF domain, but the wave functions of the 
higher band of the valence bands have large amplitude 
on the inner sites of AF domain. So the lower band has 
a large band width because of the transfer energy, and at 
M-point(fc = (7r/2,7r/2)) its energy become higher. But 
at V = 1.5, the wave function which has large amplitudes 
on the inner sites have larger energy because of the direct 
nearest-neighbor repulsion. Therefore for larger V , the 
energy difference of these band becomes large. Also the 
same reason that the midgap state becomes flat, causes 
the band invariant under the variation of V. 
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Fig. 15. Energy band of the diagonal stripe state at in T = 0.1 
and V = 0.0(a) and 1.5(b). The inset is the first Brillouin zone. 
fj, is the chemical potential. 



We show the energy band of the CDW-AF state at 
T = 0.5 and V = 1.0 and 3.0 in Fig. 16. There are two 



kinds of the bands, almost occupied valence bands and 
empty conduction band. Both kinds of the bands have 
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Fig. 16. Energy band of the CDW-AF state at in T = 0.5 and 
V = 1.0(a) and 3.0(b). The inset is the first Brillouin zone, fi is 
the chemical potential. 



two separated band. At V = 1.0 the wave function of 
the higher energy one of valence bands has large ampli- 
tude on the site where large spin exists, with the same 
spin. And that of the lower energy one has large ampli- 
tude on the site where a small spin exists, with the same 
spin. Therefore the energy difference comes from the 
direct nearest neighbor Coulomb interaction. For con- 
duction bands, also the energy difference comes from the 
direct nearest neighbor Coulomb interaction. The band 
gap between conduction bands and valence bands comes 
from the on-site Coulomb interaction. For V = 3.0, the 
split of two of valence bands becomes large. The wave 
functions of the lower band of the valence bands have 
large amplitude on the site where large spin exists, with 
the same spin, though that of the upper band have large 
amplitude on the same site but with the opposite spin. 
Therefore the energy difference of these two bands comes 
from the on site Coulomb interaction. The wave func- 
tions of both of the conduction bands have large ampli- 
tude at the site where the charge density almost van- 
ishes. Therefore the band gap energy comes from the 
nearest neighbor Coulomb interaction, and the on-site 
Coulomb interaction does not affect much on the ener- 
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gies of these two conduction bands. As we have shown, 
the band structure of the CDW-AF state varies strongly 
with the nearest neighbor Coulomb interaction. 

We show the energy band of the vertical stripe CDW- 
AF state at T = 0.1 and V = 1.75 and 2, 75 in Fig. 17. 
At V = 1.75, there are midgap bands with valence bands 




Fig. 17. Energy band of the vertical CDW-AF state at in T = 0.1 
and a)V = 1.75 and b)V = 2.75. The inset is the first Brillouin 
zone, fi is the chemical potential. 



and conduction bands. The wave function of valence 
bands has large amplitude with the same spin as the spin 
density of the site, independent of the magnitude of the 
spin density. The wave function of the conduction bands 
is opposite. But for midgap bands, the wave fuction 
has large amplitude at the site where the spin density 
is small with the same spin. The band gap energy of 
valence bands and conduction bands comes from the on- 
site Coulomb interaction, similar to the CDW-AF state 
at V = 1.0. 

But at V — 2.75, the midgap bands almost merge to 
valence bands and this phase becomes metallic unlike to 
other stripe phases. Similar to the CDW-AF states, the 
wave function of the valence bands has large amplitude at 
the site where the charge density is large, and that of the 
conduction band has large amplitude at the site where 
the charge density is small. So, the band gap energy 
comes from the nearest neighbor Coulomb interaction. 
This vertical stripe CDW-AF state changes largely the 



band structure when the value of V changes just 1. This 
feature is different from the diagonal stripe state of which 
band structure is invariant to the variation of V. 

§4. Conclusion 

In this paper, we have solved self-consistent equation 
of the extended Hubbard model on the square lattice. We 
have obtained the nearest neighbor Coulomb repulsion 
V vs. the temperature T phase diagram. At low tem- 
perature there are stripe phases and the diagonal stripe 
state is stable for the finite but small V. For larger V, 
the vertical stripe CDW-AF state becomes stable. The 
band structure of the diagonal stripe state is invariant 
against V, but the band structure of the vertical stripe 
CDW-AF state changes largely by small variation of V. 

In relation to the real materials, especially nickelates, 
the stripe charge ordered state is not stable in this 
parameter region within the mean-field approximation. 
The nearest neighbor Coulomb interaction V only sta- 
bilizes the CDW based states. In order to obtain the 
stripe charge ordered state, we think that we must take 
account the spin fluctuation, because the mean-field ap- 
proximation overestimates the magnetic order. 
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